LANDAU'S HALF-SPACE PROBLEM OF DEGENERATE 
PLASMA OSCILLATIONS WITH SPECULAR - 
ACCOMMODATIVE BOUNDARY CONDITIONS 



4 



N.V. GritsienkotI, A. V. Latyshe\o and A. A. YushkanoMj 



Faculty of Physics and Mathematics, 
Moscow State Regional University, 105005, 
Moscow, Radio st., 10- A 



In the present paper the hnearized problem of half-space plasma oscillations 
in external longitudinal alternating electric field is solved analytically. Specular ~ 
accommodative boundary conditions of electron reflection from the plasma bound- 
ary are considered. Coefficients of continuous and discrete spectra of the problem 
are found, and electron distribution function on the plasma boundary and electric 
field are expressed in explicit form. 

Refs. 34. Figs. 2. 

Keywords: degenerate plasma, half-space, normal momentum accommodation co- 
efficient, specular - accommodative boundary condition, plasma mode, expansion 
by eigen functions, singular integral equation. 

PACS numbers: 52.35.-g, 52.90.+Z 

1. INTRODUCTION 

The present paper is devoted to degenerate electron plasma behaviour 
research. Analysis of processes taking place in plasma under effect of ex- 
ternal electric field, plasma waves oscillations with various types of condi- 
tions of electron reflection from the boundary has important significance 
today in connection with problems of such intensively developing fields as 
microelectronics and nanotechnologies |[I] - 
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The concept of " plasma" appeared in the works of Tonks and Langmuir 
for the first time (see |7]-P]), the concept of "plasma frequency" was 
introduced in the same works and first questions of plasma oscillations 
were considered there. However, in these works equation for the electric 
field was considered separately from the kinetic equation. 

A. A. Vlasov |ID] for the first time introduced the concept of "self- 
consistent electric field" and added the corresponding item to the kinetic 
equation. Now equations describing plasma behaviour consist of anchor 
system of equations of Maxwell and Boltzmann. The problem of electron 
plasma oscillations was considered by A. A. Vlasov lIDj by means of solution 
of the kinetic equation which included self-consistent electric field. 

L.D. Landau had supposed that outside of the half-space contain- 
ing degenerate plasma external electromagnetic field causing oscillations 
in plasma is situated. By this Landau has formulated a boundary condi- 
tion on the plasma boundary. After that the problem of plasma oscillation 
turns out to be formulated correctly as a boundary problem of mathemat- 
ical physics. 

In jTT] L.D. Landau has solved analytically by Fourier series the prob- 
lem of collisionless plasma behaviour in a half-space, situated in external 
longitudinal (perpendicular to the surface) electric field, in conditions of 
specular reflection of electrons from the boundary. 

Further the problem of electron plasma oscillations was considered by 
many authors. Full analytical solution of the problem is given in the works 
and . 



This problem has important significance in the theory of plasma (see, 
for instance, [2], [H] and the references in these works, and also [IB], [IS])- 

The problem of plasma oscillations with diffuse boundary condition was 
considered in the works [E], [IB] by method of integral transformations. In 
the works [IS], [2D] general asymptotic analysis of electric field behaviour 
at the large distance from the surface was carried out. In the work [12] 
particular significance of plasma behaviour analysis close to plasma reso- 
nance was shown. And in the same work [ID] it was stated that plasma 
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behaviour in this case for conditions of specular and diffuse electron scat- 
tering on the surface differs substantially. 

In the works [21] and [22] general questions of this problem solvability 



were considered, but diffuse boundary conditions were taken into account. 
In the work [22] structure of discrete spectrum in dependence of parameters 



of the problem was analyzed. The detailed analysis of the solution in 
general case in the works mentioned above hasn't been carried out conside- 
ring the complex character of this solution. 

The present work is a continuation of electron plasma behaviour in 
external longitudinal alternating electric field research [i21j - [26] . 

In the present paper the linearized problem of half-space plasma oscil- 
lations in external alternating electric field is solved analytically. Specu- 
lar - accommodative boundary conditions for electron reflection from the 
boundary are considered. In [21]-[2S] diffuse boundary conditions were 
considered. 

The coefficients of continuous and discrete spectra of the problem are 
obtained in the present work, which allows us to derive expressions for 
electron distribution function at the boundary of conductive medium and 
electric field in explicit form, to reveal the dependence of this expressions 
on normal momentum accommodation coefficient and to show that in the 
case when normal electron momentum accommodation coefficient equals 
to zero electron distribution function and electric field are expressed by 
known formulas obtained earlier in [12], [T3] . 

The present work is a continuation of our work [27], in which questions 
of plasma waves reflection from the plane boundary bounding degenerate 
plasma were considered. 

Let us note, that questions of plasma oscillations are also considered in 
nonlinear statement (see, for instance, the work [2B], [2S])- 

2. PROBLEM STATEMENT 



Let degenerate plasma occupy a half-space a: > 0. 

We take system of equations describing plasma behaviour. As a kinetic 
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equation we take Boltzmann — Vlasov r-model kinetic equation: 

f + v|^ + eE|^ = M^iW(£:M. (1.1) 
ot or op T 

Here / = /(r, v, t) is the electron distribution function, e is the electron 
charge, p = mv is the momentum of an electron, m is the electron mass, 
r is the character time between two collisions, E = E(r,t) is the self- 
consistent electric field inside plasma, f^q = /eg(r, t) is the local equilibrium 
Fermi — Dirac distribution function, f^q = 0(£i?(t, x) — 8,), where Q{x) is 
the function of Heaviside, 



e{x) = 



2.p{t,x) = ^'mvp{t,x) is the disturbed kinetic energy of Fermi, £ = ^mv'^ 
is the kinetic energy of electron. 

Let us take the Maxwell equation for electric field 

divE(r,t) =47r/)(r,t). (1.2) 

Here p{r, t) is the charge density, 

p(r, t) = ej (/(r, v, t) - /o(v)) d^F, (1.3) 

where 

2d p 

= TTTT^v = dpxdpydp^. 

Here /o is the undisturbed Fermi — Dirac electron distribution function, 

/o(£) = e(£i.-£), 

h is the Planck's constant, v is the effective frequency of electron collisions, 
z/=l/r, £i^ = ^mvp is the undisturbed kinetic energy of Fermi, vp is the 
electron velocity at the Fermi surface, which is supposed to be spherical. 

We assume that external electric field outside the plasma is perpendi- 
cular to the plasma boundary and changes according to the following law: 
Eq exp{—iujt) . 
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Then one can consider that self-consistent electric field E(r, t) inside 
plasma has one ic-component and changes only lengthwise the axis x: 

E = {E^{x,t),0,0}. 

Under this configuration the electric field is perpendicular to the boun- 
dary of plasma, which is situated in the plane x = 0. 

We will linearize the local equilibrium Fermi — Dirac distribution feq 
in regard to the undisturbed distribution /o(S): 

/e, = /o(£) + Mx,i)-£]^(£F-£), 

where 6{x) is the delta - function of Dirac. 

We also linearize the electron distribution function / in terms of absolute 
Fermi — Dirac distribution /o(£): 

/ = /o(£) + /i(x,v,t). (1.4) 

After the linearization of the equations (1.1)-(1.3) with the help of (1.4) 
we obtain the following system of equations: 

dfi ^ dfi 

^ + ^.^ + ./i(^,v,t) = 

= diEp - £) [eE,(x, t)v, + iy[e.Fix, t) - 8^.]] , (1.5) 
dEg^i^.t) _ Sire f > . ,3 / 

From the law of preservation of number of particles 

J feqdnp = J fd^F 

we find: 

[8,F{x,t) - Ef] J 6{8,F - 8,)d'p = J hd'p. (1.7) 

From the equation (1.5) it is seen that we should search for the function 
fi in the form proportional to the delta - function: 

h = 8,FSi8,F-8,)H{x,fi,t), fi = ^. (1.8) 

V 
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The system of equations (1.5) and (1.6) with the help of (1.7) and (1.8) 
can be transformed to the following form: 

dH dH 



Ex{x,t) + - J H{x,ii,t)dii, 



Further we introduce dimensionless function 

UC p 

and pass to dimensionless coordinate x\ = xjl ^ where / = i^i^r is the mean 
free path of electrons, and we introduce dimensionless time i\ = vt. We 
obtain the following system of equations: 

dH dH 

-wr + ^ iyH[xi,ii,ti) = 

oti ox I 

= lie{xi,ti) -\-^ J H{xi,ii ,ti)dii , (1.9) 

Here cjp is the electron (Langmuir) frequency of plasma oscillations, 

^ m 

N is the numerical density (concentration), m is the electron mass. 

We used the following well-known relation for degenerate plasma for the 
derivation of the equations (1.9) and (1.10) 



h 

Let A; to be a dimensional wave number, and let us introduce dimension- 
less wave number ki = k—, then kx = ^1^^ where e = —. We introduce 



the quantity Ui = ut = —. 

V 
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3. BOUNDARY CONDITIONS STATEMENT 

Let us outline the time variable of the functions H{xi, /i, ti) and e{xi, ti), 
assuming 

H{xi, II, ti) = e~^'^^*i/i(xi, yu), (2.1) 

e{xi,h) = e-'^^'^e{xi). (2.2) 

The system of equations (1.9) and (1.10) in this case will be transformed 
to the following form: 

1 

dh If 
fi-^ -\- {1 - iuji)h{xi, ii) = iie{xi) + - J h{xi, fi')dfi\ (2.3) 

-1 

^.g/M... (2.4) 

-1 

Further instead of xi, ti we write x, t. We rewrite the system of equations 
(2.3) and (2.4) in the form: 

1 

dh If 

^~dx ^ '^o^^^' ^ ^^^^^ ^ 2 / ^^"^^ 

-1 

1 

"^^^^^ ^ ^ hix,ii')dii'. (2.6) 



Here 



dx J 

-1 



zq = 1 — iuji = 1 = 1 — iujT. 

V 



We consider the external electric field outside the plasma is perpendi- 
cular to the plasma boundary and changes according to the following law: 
£^oexp(— ia;t). This means that for the field inside plasma on the plasma 
boundary the following condition is satisfied: 

e(0) = eo. (2.7) 
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The non-flowing condition for the particle (electric current) flow through 

the plasma boundary means that 

1 

J iih{0,i^)di^ = 0. (2.8) 

-1 

In the kinetic theory for the description of the surface properties the 
accommodation coeflicients are used often. Tangential momentum and 
energy accommodation coeflicients are the most-used. For the problem 
considered the normal electron momentum accommodation under the scat- 
tering on the surface has the most important significance. 

The normal momentum accommodation coefficient is defined by the 
following relation: 

«P = ^^, O^a^^l, (2.9) 

where Pi and Pr are the flows of normal to the surface momentum of 
incoming to the boundary and reflected from it electrons, 



Pi = J /i^/i(0, /i)d/i, (2.10) 
-1 

1 

Pr = J fi^h{0,ii)dii, (2.11) 



quantity Ps is the normal momentum flow for electrons reflected from the 

surface which are in thermodynamic equilibrium with the wall, 

1 

Ps = J fi^hsii^)dfi, (2.12) 



where the function 

hsifi) =As, < // < 1, 

is the equilibrium distribution function of the corresponding electrons. This 
function is to satisfy the condition similar to the non-flowing condition: 

1 

J iih{0,ii)dii + J iihs{^i)d^i = 0. (2.13) 
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We are going to consider the relation between the normal momentum 
accommodation coefficient ap and the diffuseness coefficient q for the case 
of specular and diffuse boundary conditions which are written in the follo- 
wing form: 

/i(0,yw) = (1 -g)/i(0,-yw) + as, 0<yw<l. (2.14) 

Here q is the diffuseness coefficient (0 ^ g ^ 1), is the quantity 
determined from the non-flowing condition. 
From the non-flowing condition we derive 

1 1 

J iihiQ^ ii)dii = j iih(0, ii)dii -\- J /i/i(0, = 0. 

-1 -1 

In the second integral we replace the integrand according to the right- 
hand side of the specular-diffuse boundary condition (2.14). After that, 
using the obvious change of integration variable, we obtain that 



as = —2q J /i/i(0, iJ,)dii. 
-1 

Let us use the boundary condition (2.13). Using the analogous to the 
preceded line of reasoning we get: 





As = -2 J iih{0, iijdji. 



-1 

From the two last equations we flnd that 



as = qAs. (2.15) 
Further we flnd the difference between two flows 

1 

Pi — Pr = J i2^h(0, iJ,)diJ, — J /i^/i(0, /i)(i/i. 

-1 

In the second integral we use the boundary condition (2.14) again. With 
the help of (2.15) we obtain that 

1 

Pi — Pr = q j iJ.^h{0, /i)d/i — J ji^agdji = 
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Substituting the expressions obtained to the definition of the normal mo- 
mentum accommodation coefficient, we have: 



. Thus, for specular - diffuse boundary conditions normal momentum 
accommodation coefficient ap coincides with the diffuseness coefficient q. 

Equally with the specular - diffuse boundary conditions another vari- 
ants of boundary conditions are used in kinetic theory as well. 

In particular, accommodation boundary conditions are used widely. 
They are divided into two types: diffuse - accommodative and specular 
- accommodative boundary conditions (see [150]). 

We consider specular - accommodative boundary conditions. For the 
function h this conditions will be written in the following form: 



If in (2.16) we assume ^40 = ^1 = 0, then specular - accommodative 
boundary conditions pass into pure specular boundary conditions. 

Coefficients Ai and A2 can be derived from the non-ffowing condition 
and the definition of the normal electron momentum accommodation coef- 
ficient. 

The problem statement is completed. Now the problem consists in find- 
ing of such solution of the system of equations (2.5) and (2.6), which sat- 
isfies the boundary conditions (2.7)-(2.13). Further, with the use of the 
solution of the problem, it is required to built the profiles of the distribu- 
tion function of the electrons moving to the plasma surface, and profile of 
the electric field. 



P^^P^ ^ qPi-qP^ 
P - P P - P 

^ I ^ S ^ I ^ s 



= q- 



h{0, 11) = h{0,-ii) + Ao + Aiii, < /i < 1. 



(2.16) 



4. THE RELATION BETWEEN FLOWS AND BOUNDARY 

CONDITIONS 
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First of all let us find expression which relates the constants Aq, Ai from 
the boundary condition (2.16). To carry this out we will use the condition 

of non-flowing (2.12) of the particle flow through the plasma boundary, 
which we will write as a sum of two flows: 

1 

No = J iJ,h{0, /ijd/j. + J iih{0, iJ,)diJ, = 0. 

-1 

After evident substitution of the variable in the second integral we ob- 
tain: 



1 



/i(0, /i) - /i(0, -n) 



d/i = 0. 



No = J II 



Taking into account the relation (2.16), we obtain that 

With the help of this relation we can rewrite the condition (2.16) in the 
following form: 

/i(0,//) = /i(0, -fi) + Aiiii - ^), < // < 1. (3.1) 

We consider the momentum flow of the electrons which are moving to 
the boundary. According to (3.1) we have: 

Pr = Pr- ^Ai. (3.2) 

It is easy to see further that 

p. = Y- (3-3) 

With the help of the formulas (3.2) and (3.3) we will rewrite the definition 
of the accommodation coefficient (2.9) in the form: 

c^pPr - c^pY + - «p) = 0- (3-4) 
Let us consider the condition (2.13). We rewrite it in the following form: 



A_. 



'■s 



+ J /j.h{0, /j.)d/j. = 0. 
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From this condition we obtain 





As = —2 J iih{0, ii)dii = 2 J iih{0, —ii)dii. 



-1 

Using the condition (3.1), we then get 

1 

'^s = 2 J idh{0, id)did. (3.5) 



Now with the help of the second equality from (2.16) and (3.5) we rewrite 
the relation (3.4) in the integral form: 

1 



j (a^^ - ^A^)^(0,/i)c?/i = --^(1 - ap)^i. (3.6) 







Now the boundary problem consists of the equations (2.5) and (2.6) and 
boundary conditions (2.7), (3.1) and (3.6). 

5. SEPARATION OF VARIABLES AND CHARACTERISTIC 

SYSTEM 

Application of the general Fourier method of the separation of variables 
in several steps results in the following substitution 



hrj{x, ii) = exp( )$(r/, /i), (4.1) 

e^{x) = eM-—)m. (4.2) 

where t] is the spectrum parameter or the parameter of separation, which 
is complex in general. 

We substitute the equalities (4.1) and (4.2) into the equations (2.5) and 
(2.6). We obtain the following characteristic system of equations: 

1 

zo(r/ - ^)$(r/, /i) = wE{ri) + | ^ $(r/, /i') V, (4-3) 
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m = -4--l I ^iv^f^W- (4.4) 



-1 



Let us introduce the designations: 



^1 ^1 2 ^^^0 

7= 1, 7>-l, 7?1 = 



.1 + 7 
zo = l- I . 

e 

Substituting the integral from the equation (4.4) into (4.3), we come to 
the following system of equations: 

ir]-fiMr],fi) = ^{w-ril), (4.5) 

Zq 

1 

-r]lEir^) = l j Hl.t^W- (4.6) 



Here 



2 .^(1 + 7) 

^1 = V - ^ 



3 3 

Solution of the system (4.5) and (4.6) depends essentially on the condi- 
tion whether the spectrum parameter rj belongs to the interval — 1 < 77 < 1. 
In connection with this the interval — 1 < 77 < 1 we will call as continuous 
spectrum of the characteristic system. 

Let the parameter 77 G (—1,1). Then from the equations (4.5) and (4.6) 
in the class of general functions we will find eigenfunction corresponding 
to the continuous spectrum: 

$(r7,/.) = F(r7,/.)^, (4.7) 

zq 

where 

f(,,M) = F^^^-c^^(.-M). (4.8) 

?7 — /i 77 

In the equation (4.8) 6{x) is the delta-function of Dirac, the symbol 
Px~^ means the principal value of the integral under integrating of the 
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expression x ^, the function A(^) is called as dispersion function of the 
problem, 

1 

A(z) = 1 + - / ^^^^rf/i. (4.9) 
c ] II - z 

-1 

Functions (4.8) are called eigenfunction of the continuous spectrum, 
since the spectrum parameter 77 fills out the continuum ( — 1, +1) compactly. 
The eigensolutions of the given problem can be found from the equalities 
(4.7). 

The dispersion function X{z) we express in the terms of the Case dis- 
persion function [ jHT] : 

Xiz) = !-- + -( l-^)x,{z), 

Zo Zo\ T]fJ 

where 



Xc(z) = 1 + 



1 1 
z f dr I f T dr 



2 J r- z 2 J r- z 

-1 -1 



is the Case dispersion function [31] . 

The boundary values of the dispersion function from above and below 
the cut (interval (—1, 1)) we define in the following way: 

^^(/^) = lim lie (-1,1). 

The boundary values of the dispersion function from above and below 
the cut are calculated according to the Sokhotzky formulas 

2rilzo 

from where 



A±(/i) = X{ii) ± ■TTYT^'^l - ^^)' -1 < Ai < 1, 



A+(/i) - A"(/i) = ^^ii{ril - li^), 
A+(/i) + A-(/i) 

— ^ = A(/i), -l<ii<l., 

where 

1 ^ ^ 

Kf^) = 1 + / — — dr], 

^mzQ J T]- n 

-1 
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and the integral in this equahty is understood as singular in terms of the 
principal value by Cauchy. Besides that, the function X(ii) can be repre- 
sented in the following form: 

6. EIGENFUNCTIONS OF THE DISCRETE SPECTRUM 

According to the definition, the discrete spectrum of the characteristic 
equation is a set of zeroes of the dispersion equation 

^ = 0. (5.1) 

We start to search zeroes of the equation (5.1). Let us take Laurent 
series of the dispersion function: 

A(^) = Aoo + ^ + ^ + --- , \z\>l. (5.2) 



Here 



Z2 z4 



Aoo - A(oo) = 1 - i + -i^. 



zq\3 br]f 
1 /I 1 



zoVS 7r]l 

We express these parameters through the parameters 7 and e: 

2j + i£ + 7(7 + ie) 



Aoo = A(cxd) = 



(l + 7 + i£)2 
9 + 5ie{l + 7 + ie) 



15(l + 7 + i£)2 ' 

_ 15 + 7^£(1 + 7 + ze) 
35(l + 7 + z£)2 • 

It is easy seen that the dispersion function (4.9) in coUisional plasma 

(i.e. when £ > 0) in the infinity has the value which doesn't equal to zero: 

Aoo = A(oo) ^ 0. 
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Hence, the dispersion equation has infinity as a zero r/j = oo, to which 
the discrete eigensolutions of the given system correspond: 

This solution is naturally called as mode of Drude. It describes the 
volume conductivity of metal, considered by Drude (see, for example, |32])- 

Let us consider the question of the plasma mode existence in details. We 
find finite complex zeroes of the dispersion function. We use the principle 
of argument. We take the contour (see Fig. 1) 



which is passed in the positive direction and which bounds the biconnected 
domain Dji. This contour consists of the circumference 

r^ = {2: \z\ = R, R=-, e > 0}, 

£ 

and the contour 7£, which includes the cut [—1,+!], and stands at the 
distance of £ from it. 

Let us note that the dispersion function has not any poles in the domain 
Dji. Then according to the principle of argument the number [33] of zeroes 
N in the domain Dji equals to: 



N = — (fdlnX(z). 
2m J ^ ^ 



Considering the limit in this equality when £ ^ and taking into ac- 
count that the dispersion function is analytic in the neighbourhood of the 
infinity, we obtain that 

1 1 

N = ^ [ d lnA+(r) [ d lnA"(r) 

2m J ^ ' 2mJ ^ ' 

-1 -1 



2711 J A (/i) 



-1 
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So, we obtained that 

27riJ A-(r) 
-1 

We divide this integral into two integrals by segments [—1,0] and [0,1]. 
In the first integral by the segment [—1,0] we carry out replacement of 
variable r —r. Taking into account that A+(— r) = A"(r), we obtain 
that 



1 1 



AT = — - d\n — f4 = — / In TZT^ = - arg ^ 



1 



1 







(5.3) 



Here under symbol arg G{t) = arg | we understand the regular 

branch of the argument, fixed in zero with the condition: argG(O) = 0. 
We consider the curve 

7 = {z: z = G(r), ^ r ^ +1}, 

where 

It is obvious that 

G(0) = 1, lim G(r) = 1. 

T— >-+l 

Consequently, according to (5.3), the number of values equals to 
doubled number of turns of the curve 7 around the point of origin, i.e. 

N = 2m{G), 

where 

x(G) = Ind[o,+i]G(r) 

is the index of the function G{t). 

Thus, the number of zeroes of the dispersion function, which are situated 
in complex plane outside of the segment [—1,1] of the real axis, equals 
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to doubled index of the function G(r), calculated on the "semi-segment" 

[o,+i]. 

Let us single real and imaginary parts of the function out. At first, 
we represent the function in the form: 



G(//) = 



where 



A(//) = 1-- + -(1-4)M/^): 



and 



Ae(/.) = l + ^n^ ^ 



2 l + 

is the dispersion function of Case, calculated on the cut (i.e., in the interval 

(-i>i)). 

Taking into account that 

.u .1 + 7 2 £c .£(1 + 7) 

zq-1 = -%- = -% , r]^ = — = — -i , 

V e 3 3 3 

I ^^ 2 (1 + 7)' . ^(1 + 7) 

(2:0-1)771 = 1 — - — , 

we obtain 

"^^^^ pn,,)+iQ+ii,Y 

where 

P^i^i) = (1 + 7)' - Ao(Ai)(s2 - 3^2) ± + 
Q±(//) = £(1 + 7)(1 + Ao(//)) ± s(/^)(£' - 3//'). 

Now we can easily single real and imaginary parts of the function G{ii) 
out: 

Here 

^(a.) = [P\ij)? + [Q+(a.)]^ = [(1 + 1? + Ao(3m' - e')~ 
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2\-i2 



-e{l + j)sY + [£(1 + 7)(1 + Ao) - s{3fi' - e')] 

giifi) = P+(/i)P-(/i) + g+(/i)g-(/i) = [(1 + 7)' + Ao(3/i2- 
-s')]' - 6\1 + ^ns' - (1 + Ao)2] - (3/.^ - 

g^iii) = P+(/i)g-(/i) - P-(/i)g+(/i) = 2s[(l + 7)'(3/i' - e')+ 

We consider (see Fig. 2) the curve L, which is defined in implicit form 
by the following parametric equations: 

L = {il,e): gi{ii;-f,e) = 0, g2in;-f,e) = 0, < < 1}, 

and which lays in the plane of the parameters of the problem (7,£), and 
when passing through this curve the index of the function at the 

positive "semi-segment" [0, 1] changes stepwise. 
From the equation ^2 = we find: 

(1 + ^)^ = _M^)MZ£!). (5.4) 
Now from the equation pi = with the help of (5.4) we find that 



e = ^/L^), (5.5) 
where „ „ 



Substituting (5.5) into (5.4), we obtain: 

7 = -1 + (5.6) 

where ^ ^ 

^^(^^ _ . + Ao(/i)(l + XM)f 



\M[s\^i) + {I + XMY] • 

Functions (5.5) and (5.6) determine the curve L which is the border if 
the domain (we designate the external area to the domain as D~) in 
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^ 2 - 



1 - 




Fig. 2. 



explicit parametrical form (see Fig. 2). As in the work [31] we can prove 
that if (7,e) G D^, then >t:{G) = Ind[o,+i]G'(/i) = 1 (the curve L encircles 
the point of origin once), and if (7, e) G then x(G) = Ind[o^+i]G(;U-) = 
(the curve L doesn't encircle the point of origin). 

We note, that in the work [3^ the method of analysis of boundary 
regime when (7, s) E L was developed. 

From the expression (3.2) one can see that the number of zeroes of the 
dispersion function equals to two if (7,^) G D^, and equals to zero if 

Since the dispersion function is even its zeroes differ from each other by 
sign. We designate these zeroes as following ±770, by ijo we take the zero 
which satisfies the condition Re r/o > 0. The following solution corresponds 
to the zero ijo 



/i^o(x, /i) = exp 
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e^„{x) = exp 




E2. 



(5.8) 



From the equalities (5.7) and (5.8) it follows that the following solution 
corresponds to zero —770: 



This solution is naturally called as mode of Debay (this is plasma mode). 
In the case of low frequencies it describes well-known screening of Debay 
[3]. The external field penetrates into plasma on the depth of r/j, 
is the raduis of Debay. When the external field frequencies are close to 
Langmuir frequencies, the mode of Debay describes plasma oscillations 
(see, for instance, pi [52|). 



7. EXPANSION IN THE TERMS OF EIGEN FUNCTIONS 



We will seek for the solution of the system of equations (2.5) and (2.6) 
with boundary conditions (3.1), (3.6) and (2.7) in the form of linear combi- 
nation of discrete eigen solutions of the characteristic system and integral 
taken over continuous spectrum of the system. Let us prove that the 
following theorem is true. 

Theorem 7.1. System of equations (2.5) and (2.6) with boundary con- 
ditions (3.1), (3.6) and (2.7) has a unique solution, which can he presented 
as an expansion by eigen functions of the characteristic system: 




(5.9) 




(5.10) 




1 




(6.1) 







1 



e{x) = E00 + Eoexp 



( 



x\ 

- + 
7?/ 



exp 



zo-)E{r]) dj]. 



(6.2) 
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Here Eq and E^o is unknown coefficients corresponding to the discrete 
spectrum (Eq is the amplitude of Debay, Ei is the amplitude of Drude), 
E{r]) is unknown function, which is called as coefficient of discrete spect- 
rum. 

When (7, £) G D~ in expansions (6.1) and (6.2) we should take Eq = 0. 
Then from (6.1) and (6.2) we obtain: 

1 

E If/ x\ 

hix,ii) = —^li^ / exp ( - zo-)F(t], iJ.)E(r]) df], (6.1') 

Zo zqJ \ 77/ 


1 

e(x) = ^00 + y exp ( - E{r]) dr]. (6.2') 


Further we will consider the following case (7,^) G D~^. 

Our purpose is to find the coefficient of the continuous spectrum, coeffi- 
cients of the discrete spectrum and to built expressions for electron distri- 
bution function at the plasma surface and electric field. 

Proof. We substitute the expression (6.1) into the boundary condition 
(3.1). We obtain the following equation in the interval < /i < 1: 

£"00 Eq 

2 II ^ ^o{l^)+ 

zq Zq 



Zo J 





F{ri, li) - F{rj, -fi)] E{rj)drj = Ai[fi-^), (6.3) 



where the following designation is introduced 

Extending the function E{r]) into the interval (—1,0) evenly, so that 
E{r]) = E{—r]), and extending the equation into the interval (—1; 1) un- 
evenly, we transform the equation (6.3) to the form: 

1 

2% + + - [ Fin, f,)E{n)dri 

Zo Zo Zo J 

-1 
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= ^1 - ^ sign f?j . (6.4) 

We substitute the eigen functions of the continuous spectrum into the 
equation (6.4). We obtain the singular integral equation with the Cauchy 
kernel in the interval (—1, 1): 

1 

/2 
'2^^^E{r])dr]- 

-1 

-2r]lzo^^^E{ii) - zqAih = -\z^Ax sign (6.5) 
8. SOLUTION OF THE SINGULAR EQUATION 

We introduce the auxiliary function 

1 

M(^) = / ^^^^Ein)dri, (7.1) 
-1 

the boundary values of which on the real axis above and below it are related 
by the Sohotzky formulas: 

M+(/i) - M-{y) = 2m{ii^ - r]l)E{ii). (7.2) 
M1M±^ . Mi,), -1 < , < +1. (7.3) 

where 

1 

J V- 1^ 

-1 

and the singular integral in this equality is understood as singular in the 
sense of principal value of Cauchy. 

With the help of the equalities (7.2) and (7.3) and the similar equalities 
for the dispersion and auxiliary functions we reduce the equation (6.5) 
to the boundary condition of the problem of determination of analytic 
function by its jump on the contour: 
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-A {fi)[M (/i) + 2Eoo/i + Eoipoifi) - zqAiij] = 
■ZQAifi[r]i — fi ) sign fi, — 1 < /i < 1. 



This equation has general solution (see [33]): 

X{z)[M{z) + Eoifiz) + 2EooZo - zqA^z] = 
1 

2 zqAi f - r/2)sign /i 

3 2r/f 20 J fi- z 

-1 

where Ci is an arbitrary constant. 
Let us introduce auxiliary function 

1 



rpf N 1 / - Til) Sign /i ^ 

2r]fzo J II- z 



-1 



Then the general solution will be written in the form: 

X{z)[M{z) + Eoifoiz) + 2^002^0 - zoAiz] 

= '^zoAiT{z) + Ciz, 
From the general solution we can easy find the function M{z): 

M{z) = -2E^z - Eoipo{z) + zqAiz 
2 , T{z) Ciz 

Let us eliminate the pole of the solution (7.4) in the infinity. We get 
that 

Ci = {2Eoo — zoAijXoo. 

Now we eliminate the poles of the solution (7.4) in the points ±779. We 
single out items which contain polar singularity in the points z = ±770 in 
the right-hand side of the solution (7.4): 

t]qz - Tjl Ciz + IzqAiT{z) 



E, 



0- 



T]Q- Z \{z] 



26 



and 



Vqz + r]l _ Ciz + IzqAiT^z) 



rjo + z X{z) 
From the last two relations one can see that poles in the points z = ±770 
can be eliminated by one equality since the functions which compose the 
general solution are uneven: 



ZqAi = 



(7.5) 



(2/3)T(77o) - XooVo 

We find the coefficient of the continuous spectra from the Sohotzky 
formula (7.2): 

1 



M+{fj,)-M-{fj,) 



(7.6) 



We substitute the expansion (6.1) for the function h{x,iJ.) into the inte 
gral boundary condition (3.6). We derive the following equation: 

1 

^^00 + Eom{T]o) + J m{T])E{T])dT] 

1 . l-oip 

= -^^0^1 • 

00 ap 

We introduced the following designations in (7.7): 



(7.7) 



and 



m{±T]o) = J (^1? -'^fj^F{±T]o,iJ,)diJ,, 


1 

m{ri) = j {^i? -^iJ^F{T],ii)dii. 



The difference between boundary values of the auxiliary function M[z) 
from the expression (7.6) can be found with the help of the general solution 
(7.4): 



E{r^) = 



'-zqA 



2m{r]^-r]l) 3 " V\+{ri) \-{ri) 



+ 
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+(2Eoo - zoAi)Xoor] 



1 



1 



(7.8) 



■A+(r/) X-{rj). 

Let us note, that when passing the positive part of the cut (0, 1) the 
functions T{z) and X{z) make jumps, which differs from each other only by 
sign. Indeed, let us represent the formula for T{z) in the following form: 

1 r 1 

1 1 

(i/i = 

ji — z fi-t z 



1 



1 



z 



11^ — 



d/i, 



or 



Tiz) = 



2t]{zq 



dfi = — 5 — 

mzQ 



jl^ — z^ 



d/i. 



-1 <r] < +1. 



This integral can be easily calculated in explicit form. On the cut it is 

calculated by the following formula: 

Now from the Sokhotzky formula for difference of boundary values we 
obtain that when < 77 < 1: 

i7Tr]{r]l — rf) 



A+(ry) - \-{ri) = A(ry) ± 

T+(77)-r-(77) = r(/7)± 

Now one can find that 



2r]1zo 
i'nr]{rf' — r]l) 



2r]lzQ 



mr]{rf - rjl) 



2r]lzQ 



X-{tj) - X+{rj) = 2- 



2rilzQ 



Let us introduce an integral: 



Uz) = 



2r]lzQ 



9 2 

T]^ — !]( 
T] — Z 



drj. 
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It is obvious, that in the comlex plane, this integral is calculated by the 
formula: 



1 

- + ^ + (^^ 



.?)ln(i-l) 



With the help of this function we present the dispersion function in the 
following form: 

Xiz) = l-zTo{z) + zTo{-z), 
we also express the function T{z) in terms of this integral: 

T{z)=zTo{z) + zToi-z). 

Sum of the two last expressions equals to: 

X{z) + T{z) = l + 2zTo{-z). 

We note, that the integral T(—z) is not singular on the cut < rj < 1. 
The sum A(?7) + T{r]) on the cut < 77 < 1 is calculated in explicit form 
without applying to integrals: 



Kl) + TiTj) = 1 + 



1 



2r]fzo 



7] 



From the relations (7.6) and (7.8) we find the coefficient of continuous 
spectra in explicit form: 



or 



1 



2rilzoX+iri)\-ir]) 



2Eoo\oorf+ 



+^0^1377(1 + 2zTo{-z)) - X^rf 



1 



^^"^^ 277f^oA+(77)A-(77) 



{2Eoo - l)Aoo77'+ 



+20^1 (^7/(1 + 22:To(-r/))) 



(7.9) 



Let us calculate the integrals m(±77o) and 771(77) explicit form. The 
integrals 777,(±77o) can be calculated easily: 



( 4 + "») + ("0 - ^») (1 - i) 
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(-J-%) + (%^ + 5%)ln(l + i) 



Let us calculate the integral m^rj). We have: 

1 

Mr]-r]f 



P— - 2r]lzo'-^Sir] - //) 



f]-IJ, 



dji = 



1 

J (m^ - ^)") ill - 11^) 



-^-2rfMri-\)Kri)0Av)- 

II-T] 3 



Here ^+(77) is the characteristic function for the interval < ?] < 1, i.e. 

I 0, -1<77<0. 

From here we get that 

mij]) = {rf - 1)1) - ^) - ijf - Vl) [r]^ - ^v) In (l + ^) + 

+2(7^2 - r]lzo) {v-l), 



and when — 1 < 77 < 



m{r]) = {rf - rjf] 



6 



T] 



2 2 \ / 1 



The last two formulas can be written as one: 



m{T]) = {rf - rji) 



1 

L6 



77 



+ 



where 



, 1 + 

In < 77 < 1, 



In 



1 — 77 
77 



-1 < 77 < 0. 



We substitute the expansion (6.2) to the condition for the field (2.7), 
then this condition can be transformed to the following form: 

1 

£^00 + ^0 + j E{ri)dr] = eo. (7.10) 
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Singular integral equation (6.5) is solved completely. The unknown 
function E{rj) is found unambiguously and is determined by the equality 
(7.9). The unknown constants can be found from the equations (7.5), (7.7) 
and (7.10). 

9. COEFFICIENTS OF THE CONTINUOUS AND 
DISCRETE SPECTRA 

Now for finding of the coefficients of continuous and discrete spectra we 
construct system of equations (7.5), (7.10), (7.9) and (7.7): 

(2/3)r(,o) - ' ^'-'^ 

1 

Eoo + Eq + j E{r])dr] = eo, (8.2) 







or 



Eoo^Eo + ^J E{f])dr] = eo, (8.2) 



-1 



'2 



+zoAi(-7y(l + 2^To(-77)))], (8.3) 
1 

^Eoo + Eom^rjo) + J m{r])E{r])dr] 



= -l^oAi^-^. (8.4) 
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From (8.1) we find the coefficient corresponding to the mode of Debay: 
_ (2Eoo - ZQAi)\oo'nQ + lzQAiT{r]Q) 



We substitute the expressions (8.3) and (8.5) into (8.2). From here we 
derive: 
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1 - 



2Aoo77o 



(2/3)T(//o) , 



+ 



+ 



A. 



where the following designations are introduced 

1 



-J^ 



Ji = 



1 



27ri 



1 



1 



-1 



A+(77) A- (77) 



9 2 ' 
rj^ — rjf 



and 



J2 = 



2m J LA+(7y) A- (77) 



drj 



^2 _ ^2 



= eo, 



The integrals Ji and J2 can be calculated analytically with the help of 
methods of contour integration. We have 



Res +Res +Res +Res +Res 

00 -r]i rji -rjo rjo 



Z 



1 1 

+ ^ + 



27^0 



Aoo ' Ai ' X>[r^,)[r^l-4Y 



J2 = 

2r(r/i) 



Res + Res + Res + Res + Res 

00 -r]i rji -rjo rjo 



Tiz 



\{z){z^-4) 



1 



+ 



2Air/i X'irjoM " ll) 2Aic X'irjoM " ll) ' 
We substitute the integral expressions received into the preceding rela- 
tion. Therefore we obtain 



E, 



A 



00 



00 



Ai 



+ 2:0^1 



1 



A. 



6cAi 2 VAi A 



1 



1 



00 



eo, 



from where 



E, 



^00 



+ 



ZqA^ 



Ai 2Ai 



^ — (Aoo - Ai) 

6C 



Noticing that 



Aoo — Ai + 



3c' 
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we find the following relation 

^00 = eo^. (8.6) 

Substituting the deduced expression (8.6) into the equality (8.5), we 
find the following relation connecting Eq and ZqAi 

oo^O 

° ^^'{m){lo-li) ° ^ A'(/7o)(77o-^i) 

or 

Eq = -2eoAia(77o) - zqAi (-6(770) - Aooa(7?o)) , (8.7) 

where 

Let us rewrite the expression for the continuous spectrum coefficient in 
the following form 

° H3cA+(77)A-(77) 2cA+(77)A-(77)J cA+(r7)A-(77)' ^ ^ ^ 

where 

T2 = l + 277To(-77). 
We substitute (8.7), (8.8) and (8.6) into (2.4). We get the equation 

- 2eoAia(77o)m(77o) - z^Ax(h){ri^ - Xooa{'no)^m{r]o)-\- 



+Z0A1 

where 



-P2 - —Pi 
L3 2 



+ eoXiPi + zoAi^—^ = 0, 
Soap 



1 1 
If r]'^m{r])dr] If r]T2{r])m{r])dr] 

'^cj A+(/7)A-(/7)' '^ci A+(/7)A-(/7) ' 


From the last equation we find unknown constant Ai. 

2a(77o)m(7]o) - ^ - Pi 



zqAi = eoAi , . ; . (8.9) 

[\ooa{m) - iK^o)]m(77o) + §^^2 - ^^i + ^ 



Thus, all the coefficients of the expansions (6.1) and (6.2) are found 
unambiguously: the coefficient £^00 is determined according to (8.6), the 
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coefficient Eq is determined according to (8.7), the coefficient of the con- 
tinuous spectrum E(r]) is determined according to (8.8), the constant Ai 
is determined from (8.9). Finding of the coefficients of discrete and con- 
tinuous spectra of the expansions (6.1) and (6.2) completes the proving of 
existence of this expansions. Uniqueness of the solution in the form of ex- 
pansions (6.1) and (6.2) can be proved easily with the use of contraposition 
method. 

10. CONCLUSION 

In the present work new boundary conditions for the questions of plasma 
oscillations in half-space degenerate plasma were proposed. These bound- 
ary conditions are naturally called as specular accommodative conditions. 
Such boundary conditions are most adequate for the problems of normal 
propagation of plasma waves (perpendicular to the boundary), since ac- 
commodation coefficient under such boundary conditions is normal elec- 
tron momentum accommodation coefficient. 

In the present paper the analytical solution of the problem of plasma os- 
cillations in half-space degenerate plasma with accommodation of electron 
normal momentum, the coefficients of the discrete and continuous spectra 
of the problem are found in explicit form. 
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